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ABSTRACT 

Eclipsing binaries offer a unique opportunity to determine lundamcntal physical parameters of stars 
using the constraints on the geometry of the systems. Here we present a reanalysis of publicly available 
two-color observations of about 6800 stars in the Large Magellanic Cloud, obtained by the MACHO 
project between 1992 and 2000 and classified as eclipsing variable stars. Of these, less than half are 
genuine eclipsing binaries. We determined new periods and classified the stars, 3031 in total, using the 
Fourier parameters of the phased light curves. The period distribution is clearly bimodal, reflecting 
refer to the separate groups of more massive blue main sequence objects and low mass red giants. The 
latter resemble contact binaries and obey a period-luminosity relation. Using evolutionary models, 
we identified foreground stars. The presented database has been cleaned of artifacts and misclassified 
variables, thus allowing searches for apsidal motion, tertiary components, pulsating stars in binary 
systems and secular variations with time-scales of several years. 

Subject headings: binaries: eclipsing — galaxies: individual (Large Magellanic Cloud) — stars: statis- 
tics — stars: variables: other 



1. INTRODUCTION 

The last decade witnessed the birth of a new research 
field, the large-scale study of variable stars in exter- 
nal galaxies. This has been made possible by the huge 
databases of microlensing obs ervations in the Ma gellanic 
Clouds, such as the M ACHO (lAlcock et al.lli993lb OGLE 
(|Udalski et al.l Il99"l and EROS (|Aubourg et all fl995h 
projects. These programs (beyond their primary pur- 
pose) resulted in the discovery of thousands of eclipsing 
binaries and many other types of variable stars and gave 
an unprecedented homogenous coverage of their light 
curves. Recent all-sky surveys also giv e a good oppor- 
tunit y to study large numbers of stars (jPaczvhski et al.1 
2006). However, processing the huge amount of data can 
be quite challenging, even when looking at seemingly 
simple issues such as classification of variables, period 
determination, etc. 

The astrophysical potential of large databases of eclips- 
ing binaries has been explored by a number of authors, 
but we are still far from a full exploitation. Catalogs of 
eclipsing variables in the Magellanic Clou ds have been 
publi s hed with an increasing complete ness (jAlcock et al.l 
119971 IWvrzvkowski et all [20031 12004( 1. Attempts have 
been made to improve our understanding of the forma- 
tion of co ntact binaries, e sp ecially those with giant com- 
pone nts (Rucinski 1997a.b, [19981 : M aceroni fc Rucins"kl 
119991 : [Rucinski & Maceroni 2001). To help analyse large 
samples of stars, v ariou s automated pipelin es have been 
devel oped (|Devorll2005t iTamuz et aiTl2006t iMazeh et all 
2006), while the large-number statistics helped investi- 
gate orbital circularization over a broad range of stellar 
parameters (|Faccioli et alj |2005). Another application is 
determining accurate distances using early-type detached 
eclipsing binaries, which has been shown to offer the 
most accurate calibration of the local distance scale (e.g., 
iGuinan et aHHoM iRibas et al J 1200(1 : iWvithe fc Wilson! 

1 E-mail: derekas@physics.usyd.edu.au 



200lHSalaris fc Groenewegenl l2002t iFitzpatrick et all 
2002; Ri bas et alJl2002HMichalska fc Pigulskill2005H . 

In contrast to these applications, the near-decade 
long time coverage of the available data has rarely 
bee n utilised. For examp le, iPalen fc Armstrong! (|2003f ) 
and iJohnson et al.l (|2004l ) reported on searches for ter- 
tiary companions of eclipsing binaries in the MACHO 
database, but no definite results have been published. 

Here we present an analysis of the publicly avail- 
able MACHO light curves of LMC variables classified 
as eclipsing binaries, which have been online since 2001. 
The main aim of this project is to measure period changes 
and search for eclipsing binaries with pulsating compo- 
nents. The period-luminosity (P-L) relations of contact 
binaries and their relation to pulsating red giants were al- 
ready discussed in lDerekas et al.l (|2006f) . In this paper we 
discuss the general properties of the sample, presenting a 
full reclassification, newly determined periods, the color- 
magnitude diagram and period-luminosity distributions. 

2. PERIOD DETERMINATION AND CLASSIFICATION 

The MACHO observations were carried out between 
1992 and 2000 with the 1.27m Great Melbourne tele- 
scope at Mount Stromlo Observatory, Australia. The 
telescope was equipped with a specificly designed cam- 
era, which gave a 0.5 square degree field of view. The ob- 
servations were obtained in two non-standard bandpasses 
simultaneously: a 440-590 nm MACH O "blue" fil ter and 
a 590-780 nm MACHO "red" filter (|Cookl [l995h . The 
typical number of observational points obtained is about 
1200. Some of the data are publicly available on the 
MACHO website 2 , where one can select samples based 
on an automatic variability type classification. Using the 
web-interface, we individually downloaded light curves 
of all stars classified as eclipsing binaries, 6833 in total. 
However, it became obvious very quickly that the auto- 
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matic classification was not perfect, as a large fraction of 
'eclipsing binaries' turned out to be Cepheids, RR Lyrae 
stars or long-period variables. We therefore needed to 
re-classify all stars. Furthermore, besides the problems 
with the classification, we also found that the MACHO 
periods were incorrect for significant number of stars. For 
this reason, we redetermined periods for all the 6833 ob- 
jects and classified them based on the light curve shapes. 
Since the MACHO light curves contain more blue points 
than red ones, we used the former for determining peri- 
ods. 

For an initial period determination, we us ed the Phase 
Dispe rsion Minimization method (PDM, Stcllingwerf 
119781) . The method is based on a technique to mini- 
mize the sample variance of the phase diagram. It is 
calculated in the following way. The phase diagram for a 
given trial period is divided into m segments, each con- 
taining nj (j = l...m) points. The variance for each seg- 
ment is defined as s 2 = Ya=i ( m i — 5Tj) 2 / (nj — 1), where 
mj is the observed magnitude, TTTj is the mean magnitude: 
THj = YaLi m i/ n j- The sample variance s 2 = Yl^Lx s f ^ s 
minimized to get the best period. We calculated 600 
000 trial phase diagrams for each star, covering a wide 
range of periods between 0.085 days and 1000 days (from 
0.001 cycles/day to 12 cycles/day, with equidistant steps 
in frequency). 

However, this was not enough for determining the true 
periods. As commonly happens in finding periods from 
light curves of eclipsing binaries, a significant fraction of 
the PDM results were harmonics or subharmonics of the 
true period (PpDM/Ptruc was a ratio of small integers). 
To correct for this, we proceeded as follows. A visual 
inspection of every phase diagram showed whether the 
actual period was an alias or just slightly inaccurate. In 
the case of an alias, we multiplied the PDM period by 
different constants (in most cases by 2) until the shape of 
the curve was consistent with that of an eclipsing binary. 

We next us e d th e Strin g- Lengt h method 
(|Lafler fc Kinmanl H96l iClarkd [2002) to im- 
prove period determination. With this method, 
one calculates the total length of the phase dia- 
gram for any given period as the following sum: 

SL 2 = J2t=i (( m i+i _ m i) 2 + (^i+i _ </>i) 2 ) > where 
{<^i,m;}" =1 is a folded dataset sorted in phase; ipi and 
mj are the phase and the magnitude of the observations 
taken at time t; (so that ip\ = [(ti — t ) /P], where to 
was chosen as an epoch of the deeper minimum, P is the 
period, [ ] is the fractional part). The best period here 
was the one that minimized SL 2 . In our case, we applied 
the SL method for 1000 periods within ±1% of the best 
PDM period. An example of the period improvement 
is shown in Fig. [1] The typical period improvements 
resulted in a change in the 5-6th decimal place. 

During the visual inspection of the phase diagrams 
we also made a rough classification of all 6833 vari- 
ables. Based on the light curve shape, phased with the 
finally adopted periods, we placed each star into one of 
the following categories: Algol-type, (3 Lyrae-type, W 
UMa-type eclipsing binary, pulsating star (including RR 
Lyraes, Cepheids, Miras, etc.), non-periodic or multiply 
periodic, and unidentifiable. In cases of a pure sine- wave, 
which can be observed in a contact binary or a pulsating 
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Fig. 1. — An example (81.9006.28) for improvement in period de- 
termination using PDM (top panel) and SL (bottom panel) meth- 
ods. 



star, we compared the amplitudes of the MACHO blue 
and red phase diagrams. If there were noticeable color 
variations exceeding a few hundredth of a magnitude, we 
considered the star as pulsating variable. 

After the whole procedure, we ended up with 3031 gen- 
uine eclipsing or ellipsoidal binaries, the rest belonging 
to other type of variables. The identification numbers, 
J2000 coordinates, calibrated colors and magnitudes, pe- 
riods and epochs of minimum light for the eclipsing vari- 
ables are listed in TableQ] (available in full electronically). 
In the last column of the table, we note cases where a 
very faint secondary eclipse might possibly exist, which 
would result in a 1/2 ambiguity in the period. 12 of 
the 3031 systems have no periods because they showed 
signs of only one eclipse over the eight years of MACHO 
observations. 

A comparison of our periods with those in the online 
MACHO catalog shows that there is a major disagree- 
ment. A plot of the ratios of the periods (Fig. [2]) re- 
veals an interesting structure. A minor fraction of the 
period ratios are scattered around 1, while the majority 
are dominated by the cases when the MACHO period 
is half of ours. Also, we identified a number of stars 
with period ratios at 1/3, 1/4, 1/5, etc. The diagonal 
line in the lower part of the histogram corresponds to 
cases when the MACHO frequency (1/Pmacho 

) is ±1 

c/d offset from the true frequency or one of its integer 
harmonics. 

To illustrate the improvement of periods in this paper 
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Fig. 2. — A comparison of periods presented in this study and 
those of the MACHO project. 
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Fig. 3. — The ratio of the length of phase diagrams calculated 
with period in this paper to that using the MACHO period versus 
period of this paper. 
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Fig. 4. — The histogram of period ratios presented in this study 
and by the MACHO project. Note the vertical axis is in logarithmic 

goi e (EA), (3 Lyrae (EB) and W UMa (EW) type stars. 
However, this kind of classification contains some subjec- 
tivity, so we d ecided to use a more objective method for 
this purpose. iRucinskil (|1993h showed that light curves 
of W UMa systems can be quantitatively described using 
only two coefficients , &2 and &4, of the cosine decomposi- 
tion a i cos(27ri(p). iPoimah ski (2002j) tested the behav- 
ior of semidetached and detached configurations in the 
&2 — a4 plane by decomposing theoretical light curves 
into Fourier coefficients. He found that in most cases 
contact, semi-detached and detached configurations can 
be distinguished quite accurately. We chose this method 
because it can be easily implemented for a large set of 
light curves. 

Following the definitions bv lPojmariskil (|2002h . we de- 
composed every light curve in the following form: 



m(ip) = mo 



(ai cos (2mip) + hi sin (2m(p)) (1) 



compared to those in the online MACHO catalog, we 
performed the following statistical evaluation. For each 
star, we calculated the string-length using both our pe- 
riod and that of the MACHO online catalog. The ratio 
of the two lengths is plotted as a function of period in 
Fig. [3] We see a considerable improvement for all but a 
few stars; the points over 1.0 all refer to datasets with 
very few points, for which none of the statistics is well 
defined. 

Ratios other than 0.5 occurred almost exclusively in 
stars with eccentric orbits, which causes a shift in the 
secondary minimum from 0.5 phase. These problems 
clearly illustrate that one has to be very careful when 
determining periods automatically for eclipsing binaries. 

The histogram of the period ratios (Fig. [4| shows that 
roughly 16% of the periods agree. The true period turned 
out to be the double of the given MACHO period in 
about 78% of the binary sample, while the remaining 
6% have other ratios. In Fig. [5] we show a few typical 
examples, plotting phase diagrams with the MACHO and 
our periods. 

With the corrected phase diagrams, it became possi- 
ble to re-classify the sample. As discussed above, pul- 
sating and non-periodic variables were easily excluded, 
while eclipsing binaries were visually pre-classified as Al- 



where m(ip) is the phased light curve, tuq is the mean 
magnitude, while the zero point of the phase corre- 
sponds to the primary minimum. The resulting distri- 
bution in the a2 — &4 plane, based on the MACHO blue 
data, is sh o wn in Fig. [6] Using the boundary lines of 
IPoimahski (|2002f ), we marked the contact (EC), semi- 
detached (ESD) and detached (ED) configurations with 
different symbols. 

The sample is dominated by bright main-sequence de- 
tached (1482 stars) and semidetached (937 stars) bina- 
ries, while a small fraction consist of contact systems 
(612 stars). Many of the latter are foreground Milky 
Way objects, as can be shown from their positions in 
the color-magnitude diagram (Section 3.2). We list the 
EC/ESD/ED classification for each star in Table [TJ 

3. GENERAL PROPERTIES 

After the visual inspection, we identified 3031 stars 
as eclipsing binaries, which is about 44% of the down- 
loaded sample. We emphasize that our sample probably 
does not contain all LMC eclipsing binaries observed by 
the MACHO project. We assume that there might be 
eclipsing binaries classified as other variable types (such 
as RR Lyraes, Cepheids. Semiregulars, etc.). For exam- 
ple, iFaccioli et al.l (|2005f ) studied orbital circularization 
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TABLE 1 

Representative lines from the full Table 1 (available electronically). The (*) symbol in the Type column denotes 

foreground objects. 



ID number 


RA 


Dec. 


V 


V-R 




Period 


Epoch 


Type Note 






(J2000.0) 


(J2000.0) 


(mag) 


(mag) 




(days) 


(2400000+) 




47. 


1402.116 


04:49:21.544 


-68:12:31.70 


18 


.644 


-0.034 


1 


.761188 


50249.7483 


ESD 


47. 


1527.178 


04:49:28.952 


-67:55:24.15 


19 


.223 


0.725 


12 


.556414 


50244.6082 


EC 


47. 


1528.41 


04:49:34.577 


-67:49:48.88 


17 


.209 


-0.179 


5 


.024733 


50249.8034 


ED 


47. 


1529.60 


04:49:44.988 


-67:48:57.28 


18 


.005 


0.574 


112 


.073172 


50252.7272 


ED 


47. 


1531.32 


04:49:51.948 


-67:41:38.13 


16 


.842 


-0.159 





.459804 


50249.9699 


EC 


47. 


1530.30 


04:49:52.060 


-67:42:09.84 


15 


.519 


-0.249 


1 


.633056 


50249.3579 


ED 


47. 


1649.69 


04:50:16.357 


-67:53:06.87 


17 


.396 


-0.153 


1 


.453116 


50250.0195 


ESD 


47. 


1645.115 


04:50:21.514 


-68:07:37.01 


18 


.048 


-0.163 


1 


.711124 


50249.3182 


ED 


47. 


1647.143 


04:50:33.106 


-67:58:44.94 


18 


.440 


-0.186 


1 


.620856 


50250.1698 


ED 


47. 


1647.285 


04:50:34.347 


-67:59:24.49 


19 


.653 


0.082 


1 


.387449 


50249.4512 


ESD 


11. 


9720.84 


05:40:08.972 


-70:15:03.29 


17 


.506 


0.953 


4 


.607095 


50248.2811 


EC(*) 



of LMC eclipsing binaries presenting a new sample of 
4576 stars; however, that sample is not available to us. 
With this caveat, we discuss the main properties of the 
sample. 

3.1. Period distribution 

The most fundamental parameter of binary stars is the 
orbital period, whose distribution can be of great aid in 
understanding formation and evolution of close binaries. 
Eclipse detection is highly influenced by the orbital pe- 
riod, because the wider the separation of components, 
the longer the orbital period gets, implying that eclipses 
will be seen in narrower range of the inclination angle. 
Hence the chances to detect eclipses are much higher 
for short orbital periods. Wel l befo re th e microlensing 
projec ts, iFarinella fc Paolicc hi (1978) and lGiuricin et all 
(1983) found a multi-modal period distribution in the 
Milky Way. Recent studies of eclipsing binaries in the 
Small and the Large Magellanic Clouds revealed an over- 
all similarity in t he period distributi o ns, peaking between 
1 and 2 days (lAlcock et all Il997t iWvrzvkowski et all 
[2001 [200l [Devorll2005D .~ 

For a direct comparison with the OGLE sample for 
the LMC and the SMC, we downloaded periods of eclips- 
ing binaries from the OGLE internet archive 3 . Our pe- 
riod distribution (Fig. [7]) is in a good agreement with 
that of the OGLE sample of LMC and SMC stars 
(|Wvrzvkowski et all [20031 12004ft : the majority of sys- 
tems have short periods, peaking between 1 and 2 days, 
and roughly 20% of stars have periods longer than 10 
days, which is consistent for both the SMC and the LMC. 

The distribution sh own in the up per panel of Fig. [7] is 
also similar to that of iDevorl (|2005D . which was based on 
the OGLE observations of the Galactic Bulge. However, 
we do not confirm his conclusion that the peak at ~ 100 
days i s due to pulsating red giant stars. While Devor 
(2005) used single filtered I-band data of the OGLE 
project, thus having no information on the color varia- 
tions, the two-color MACHO observations clearly showed 
that those long-period variables have only slight color 
changes, which can be explained by the non-uniform tem- 
perature distribution on the surface of the non-spherical 
components. 

3 http://bulgc.astro.princeton.edu/ ogle 



The lower panel of Fig. [JJ shows the histogram of all 
eclipsi ng binaries in t he Genera l Catalog of Variable 
Stars (|Kholopov et all (|1985-1988l) with all the recent up- 
dates available for download from the GCVS website 4 ). 
Compared to the upper panel of Fig.[jj the distributions 
are different in two period ranges. For shorter periods, 
there is an excess of stars in the GCVS, which are all 
short-period main-sequence binaries that have fainter ab- 
solute magnitudes than the MACHO limit in the LMC. 
On the other hand, for periods longer than 40 days, there 
is a lack of stars in the GCVS, which we interpret as 
caused by a selection effect: these systems need years of 
observations before classification and period determina- 
tion, which was hardly possible before the development 
of automatic all-sky sur vey projects like the All-Sky Au- 
tomated Survey (ASAS, Poimanski 2002). It also implies 
that there are considerable number of bright long-period 
eclipsing binaries waiting for discovery. 

To make further comparison, we took detached and 
semi-detached binaries in the Galaxy observed by the 
ASAS project and compared their period distribution 
with detached and semi-detached binaries in the SMC 
and LMC. In the upper panel of Fig. [5]we plotted the pe- 
riod distribution of orbital periods shorter than 10 days, 
which shows an overall similarity for the three galaxies. 
This diagram also shows the presence of a selection ef- 
fect in the data: there is a noticeable dip at P=2d in the 
histograms of the MACHO and the ASAS data, which 
reveals that stars with exactly 1-d periods (i.e. the half of 
the true ones) were discarded during the initial analyses. 
For example, based on the shape of the histogram, about 
50 MACHO eclipsing binaries with P ~ 2d are missing 
from the sample, presumably due to deliberate exclusion 
from the data. 

In the lower panel of Fig. [5] we plot the cumulative pe- 
riod distribution for the three samples. While the LMC 
and SMC samples have virtually indistinguishable distri- 
butions, the Milky Way shows an excess of stars around 
10 days. A two-sample Kolmogorov-Smirnov test showed 
that the difference between the LMC and the SMC sam- 
ples is insignificant, which confirms that the MACHO 
and OGLE selections of detached and semi-detached bi- 
naries were similarly complete. On the other hand the 

4 http:/ /www. sai.msu.su/groups/cluster/gcvs/gcvs/ 
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Fig. 5. — Examples of typical period ratios. In each case, the 
left panel shows the data phased with the MACHO period, while 
the right panel shows the final phase diagram with our period. 

probability that the ASAS stars have the same distri- 
bution is only 10 _11 %. This is a simple consequence of 
the different sampling: while MACHO and OGLE data 
represent magnitude-limited samples in absolute magni- 
tudes, covering the brightest ~ 6 mag of the CMD, the 
ASAS database contains a broad mixture of stars in the 
Milky Way with a much wider range in absolute magni- 
tude. 

It is also apparent in Fig. [5] that, although the long- 
period tails of the ASAS and MACHO histograms over- 
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Fig. 6. — Classification of LMC eclipsing binaries in the Fourier 
coefficients plane a2 — a,j. 
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Fig. 7. — Upper panel: period distribution of eclipsing binaries in 
the MACHO database. Lower panel: histogram of orbital periods 
of eclipsing binaries in the General Catalog of Variable Stars. 
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P (days) 

Fig. 8. — Upper panel: Orbital period distributions of detached 
and semi-detached eclipsing binaries in the Milky Way and the 
Magellanic Clouds. Lower panel: The normalized cummulative 
distributions. 

lap very well, the short-period peak in the LMC data 
is higher because of the different distributions of stars 
in the samples. Whereas the MACHO data cover the 
upper part of the Hertzsprung-Russell diagram with a 
relatively bright limit in absolute magnitude, the ASAS 
sample contains a broader mixture of stars in the galactic 
neighborhood, so that short-period B-type systems have 
a smaller relative contribution. 

3.2. The Color- Magnitude Diagram and 
Period-Luminosity Relations 

In order to construct the Color-Magnitude Diagram 
(CMD), we converted the observed MACHO blue and red 
magnitudes into Kron-Cousin s V an d R using the equa- 
tions derived bv lAlcock et al.l (jT999) . The resulting dia- 
gram, after translating apparent magnitudes outside the 
eclipses to absolute magnitudes using //(LMC) = 18.50, 
is shown in Fig. [9] We also sh ow stellar evolutionary 
tracks bv lCastellani et ail {2003). Most of the stars are 
moderately massive near-main sequence stars, while the 
red giant branch of evolved stars is also clearly recognis- 
able. The latter is dominated by first-asce nt red giants 
and a few Asymptotic Giant Branch stars (jAlcock et al.l 
f2000h . 

In Fig. [10] we plot the CMD for four orbital period 
ranges. These plots were used for cleaning the sample 
of the foreground objects, adopting the following simple 
considerations. First, the spread of distances of eclipsing 
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Fig. 9. — The color-magnitude diagram of LMC eclipsing binaries 
with metal- poor evolutionary track s of single stars overlaid (models 
taken from lCastellani et all l|2003l \'). 

binaries in the LMC is negligible compared to the dis- 
tance to the LMC. This means that we can determine 
the absol ute magnitude for any given object wi thin ±0.1 
mag fe.g. lNikolaev et alll2004t lLah et al.ll2005h . Second, 
we assumed that the mean position of a binary star in 
the CMD can be approximated by the location of the 
brighter component. In the case of two identical compo- 
nents, this assumption means a 0.75 mag fainter abso- 
lute magnitude, however, as we proceeded, this was still 
a useful simplification in estimating foreground contam- 
ination. 

Theoretical evolutionary tracks give the basic phys- 
ical parameters of a single star, as well as color and 
magnitude information. From the physical parame- 
ters we can calculate the minimum orbital period for 
any given mass, temperature and luminosity, apply- 
ing simple relations. Combining the Stefan-Boltzmann 
law: L/Lq = (R/R e ) 2 (T/T q ) 4 and Kepler's third law: 
a 3 /P 2 = G {Ah + M 2 ) /4 • 7T 2 , we can calculate the min- 
imum orbital period of a system for two extreme cases: 

Case 1.: two identical components, for which the mini- 
mum orbital period occurs when a — 2R and Mi + Mi = 
2M; the minimal orbital period is 



Pmin(l) = 47T 



R 3 
GM 



(2) 



Case 2.: a negligible secondary component, for which 
the minimum orbital period occurs when a = R and 
Mi + M-2 = M; the minimum orbital period is 



,(1) 



R 3 
GM 



(3) 



(M is the mass of a model point of the isochrone, while 
R is calculated from L and T). 

Using Eq. ^ and Eq. and theoretical stellar evolu- 
tionary tracks, we can determine contours in the CMD, 
on which the minimum orbital period has a given value. 




Fig. 10. — Eclipsing binaries in four period ranges and estimated locations of the minimum orbital periods. Stars upper and right of the 
dashed lines have positions that are not compatible with the LMC membership. 



This is shown in the four panels of Fig. [10] by the thick 
solid lines (Case 1) and dashed lines (Case 2). As repre- 
sentative limits, we determined the locations where the 
minimum orbital period is 0.5, 1, 10 and 100 days. 

We have selected the foreground objects as follows. As 
Fig. [10] shows, moving up and right in the CMD, the 
minimum orbital period gets longer. Therefore, if a star 
with P or b < 0.5d is located above the limiting line of 
-Pmin = 0.5d, it must be in the foreground of the LMC. 



This is shown in the upper left panel of Fig. [10] As 
expected, most of these short-period reddish eclipsing 
binaries are Galactic W UMa-type stars, many magni- 
tudes above the calculated period lines (which means the 
adopted simplifications do not affect the conclusions). 

For longer periods, only a few stars are clearly fore- 
ground objects (most notably in the upper right panel 
of Fig. [TO]) , while for the majority, positions in the CMD 
are compatible with the LMC membership. We flagged 
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Fig. 11. — Period-K magnitude relation for eclipsing binaries in 
the LMC. 

all obvious foreground stars, being those located at least 
1 mag above the dashed lines, in the last column of Table 
Q] (54 stars in total). 

Finally, we briefly examine the period-luminosity (PL) 
relation of stars in the sample, for which we found near- 
infrared K- magnitudes in the 2MASS Point-Source Cat- 
alog. In iDerekas et alj (|2006h we already discussed how 
the red giant eclipsing binaries/ellipsoidal variables can 
be fitted with a simple model using Roche-lobe geometry. 
Here we examine the correlation between the period and 
the K-magnitude for the detached and semi-detached bi- 
naries. This is shown in Fig.HU where the ED/ESD/EC 
classes are plotted with different symbols. As expected, 
detached binaries are spreaded uniformly, while longer 
period semi-detached systems may follow a loose corre- 
lation (the correlation coefficient for P > lOd is ~ 0.58), 
but it is not as tight as for long-period contact binaries 
(whose correlation coefficient is ~ 0.86). Systems further 
away from the main correlation line might be foreground 
stars left unidentified as such in the CMD analysis. 

4. SUMMARY 

In this paper we present the analysis of online light 
curves of eclipsing binaries in the Large Magellanic Cloud 
monitored by the MACHO project. We downloaded the 
data of all 6833 stars classified as eclipsing binaries from 
the Variable Star Catalog Retrieval form of the project. 

We re-determined the periods for every star and re- 
classified all stars based on their light curve shape. As 
a result, 3031 stars remained as eclipsing binary, while 
the rest of the sample were RR Lyrae stars, Cepheids or 
long-period pulsators. For the binary sample, we showed 
that roughly 16% of the periods agreed with those given 
in the catalog. For almost three-quarter of the sample, 
the catalog periods turned out to be the half of the real 
ones. 

We calculated the period histogram, which shows bi- 
modal feature on logarithmic scale. The maximum of the 



distribution is at binaries with periods between 1 and 2 
days, while roughly 20% of the sample (about 600 stars) 
have period between 10 and 200 days. We compared 
our period distribution for detached eclipsing binaries, 
with those of the SMC and the Milky Way and we found 
a good agreement in the most common periods in the 
studied galaxies. 

For a more objective classification, we used cosine de- 
composition of the light curves into Fourier coefficients, 
where contact, semi-detached and detached configura- 
tions can be distinguished in a2 — a4 plane. The sam- 
ple is dominated by bright main-sequence detached and 
semi-detached systems, while a small fraction consists of 
contact binaries. However, many of the latter are fore- 
ground Milky Way objects, which can be shown from 
their position in the color-magnitude diagram. Most of 
the stars with periods shorter than 0.5 days are galactic 
W UMa systems. For longer period systems, we found a 
few obviously foreground objects and they are flagged in 
the last column of Table [TJ 

The presented database opens up a new avenue into us- 
ing the MACHO database for studying secular variations 
over a time-scale of a decade. This includes a search for 
systems with apsidal motion (recognizable through phase 
shifts of the primary and secondary minima in the oppo- 
site direction), tertiary components (indicated by cyclic 
phase shifts of the primary and secondary minima in 
the same direction) and eclipsing binaries with pulsating 
components (revealed by periodic secondary brightness 
fluctuations). The latter would be particularly interest- 
ing if evidence for tidally induced oscillations could be 
found in a larger sample of stars. With the OGLE-III 
project still taking data, after identifying these theoret- 
ically important test objects, it would become possible 
to extend the full time-coverage of observations up to 
14-15 years. In a companion paper we will discuss the 
measured period changes for the whole sample of 3031 
eclipsing binaries. To help the interested researchers, we 
have collected the original MACHO blue and red light 
curve data in a single compressed data file that can be 
accessed as an electronic appendix to this paper. 
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